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We study the amplitude modulation of ion-acoustic wave (lAW) packets in an unmagnetized 
electron-ion plasma with two-temperature (cool and hot) electrons in the context of the Tsallis’ 
nonextensive statistics. Using the multiple-scale technique, a nonlinear Schrodinger (NLS) equation 
is derived which governs the dynamics of modulated wave packets. It is shown that in nonextensive 
plasmas, the lAW envelope is always stable for long-wavelength modes {k —>■ 0) and unstable for 
short-wavelengths with fc > 1. However, the envelope can be unstable at an intermediate scale 
of perturbations with 0 < fc < 1. Thus, the modulated I AW packets can propagate in the form 
of bright envelope solitons or rogons (at small- and medium scale perturbations) as well as dark 
envelope solitons (at large scale). The stable and unstable regions are obtained for different values 
of temperature and density ratios as well as the nonextensive parameters qc and qh for cool and hot 
electrons. It is found that the more (less) the population of superthermal cool (hot) electrons, the 
smaller is the growth rate of instability with cutoffs at smaller wave numbers of modulation. 

PACS numbers: 52.35.Fp; 52.27.Ep. 


I. INTRODUCTION 

Over the last few decades, wave dynamics in multi- 
component plasmas has been a front area of research 
for the plasma community. Various types of nonlin¬ 
ear waves namely, ion-acoustic waves (lAWs), electron- 
acoustic waves (EAWs) etc. as well as nonlinear coherent 
structures associated with them namely, solitons, enve¬ 
lope solitons, shocks etc. have also been studied in dif¬ 
ferent plasma environments. As is known, ion-acoustic 
waves are electrostatic plasma modes, where a popula¬ 
tion of inertial ions oscillate against a dominant ther- 
malized background of electrons that provide necessary 
restoring force. The phase velocity of such waves lies be¬ 
tween the electron and ion thermal velocities. Plasmas 
with two groups of electrons (e.g., cool and hot electrons) 
are very common in laboratories m as well as in space 
environments [3]. In this context, a large number of in¬ 
vestigations have been made to study the characteristics 
of lAWs in two-temperature electron plasmas [MS]- Such 
two-temperature plasmas have been observed by various 
satellite missions, e.g., Fast Auroral Snapshot (FAST) at 
the auroral region (TU), Viking Satellite |3] and THFMIS 
mission m- Furthermore, several investigations made 
by GFOTAIL and POLAR [T^j in the magnetosphere 
have also confirmed the existence of such electron popu¬ 
lations. Jones et al. [5| reported experimentally the exis¬ 
tence of ion acoustic solitary waves (lASWs) in plasmas 
with hot and cold electrons. Baluku et al. m studied 
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small-amplitude ion-acoustic solitons in a plasma with 
two-temperature kappa distributed electrons and found 
that superthermality of electrons played an important 
role for the modification of solitary structures. 

The linear properties of the lAWs have been exten¬ 
sively studied and well understood in multi-component 
plasmas in presence of hot electrons obeying Maxwellian 
and non-Maxwellian velocity distributions. When the 
nonlinear effects are concerned, the formation of ion- 
acoustic solitons may be possible due to a delicate bal¬ 
ance between the nonlinearity and the dispersion which 
also has been anticipated theoretically via the deriva¬ 
tion of Korteweg-de Vries (KdV) or Zakharov-Kuznetsov 
(ZK) equations for small-amplitude perturbations [Ill 
US), and via the Sagdeev potential formalism, accounting 
for the arbitrary amplitude excitations [16] . Existence of 
these nonlinear structures was experimentally confirmed 

m- 

The space observations and laboratory experiments 
indicate the presence of particles which obey non- 
Maxwellian velocity distributions in plasmas. The ob¬ 
servations of the voyager PLS muni have also shown 
the existence of both cool and hot electrons populations 
(non-Maxwellian) in Saturn’s Magnetosphere. One of 
such distributions as reported by Renyi [20] and pro¬ 
posed by Tsallis |2l| is the Boltzmann-Gibbs-Shannon 
(BGS) entropy in which the degree of nonextensivity of 
the plasma particles is characterized by the entropic in¬ 
dex q, i.e., the distribution function with q < 1 , com¬ 
pared with the Maxwellian one {q = 1) indicates the 
system with more superthermal particles (superexten- 
sivity), whereas the q distribution with g > 1 is suit¬ 
able for plasmas containing a large number low-speed 
particles (subextensivity). However, because of long- 
range Goulomb force between plasma particles and the 
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presence of many superthermal particles in astrophysi- 
cal and space environments a g-distribution with g < 1 is 
strongly suggested for real plasmas or superthermal plas¬ 
mas. Such g-nonextensive distribution has been widely 
used to study various nonlinear coherent structures like 
solitons, double layers or shocks etc. [22H28] . 

The propagation of wave packets in a nonlinear, dis¬ 
persive medium is subjected to the modulation of their 
wave amplitudes, i.e., a slow variation of the wave packets 
envelope due to the nonlinear self-interaction of the car¬ 
rier wave modes. The systems evolution is then governed 
through the nonlinear Schrodinger (NTS) equation and 
the associated modulational instability (MI). The lat¬ 
ter in nonlinear, dispersive media has been a well-known 
mechanism for the localization of wave energy. Further¬ 
more, such instability results into the exponential growth 
of a small plane wave perturbation which leads to the am¬ 
plification of the sidebands, and thereby breaking up the 
uniform wave into a train of oscillations. Thus, the MI of 
wave packets in plasmas acts as a precursor for the forma¬ 
tion of bright envelope solitons or highly energetic rogue 
waves (rogons), otherwise the dark envelope solitons may 
be formed. The amplitude modulation of electrostatic 
waves in plasmas have been investigated by a number of 
authors owing to their importance in space, astrophys- 
ical and laboratory plasmas. A great deal of attention 
has been given to the MI oi lAWs in non-Maxwellian 
plasmas. To mention few, Saini and Kourakis [29] inves¬ 
tigated the amplitude modulation of dust-acoustic waves 
in dusty plasmas with superthermal ions. Similar inves¬ 
tigations have also been made of lAWs in plasmas with 
two-temperature electrons m and kappa-distributed su¬ 
perthermal electrons [3l]. Furthermore, Bains et al. [25] 
studied the MI of ion-acoustic wave envelopes in a plasma 
containing non-extensive electrons. 

On the other hand, the generation of highly energetic 
rogue waves (or rogons) in plasmas has been a topic of 
important research. Such waves may also appear in other 
environments such as oceanic rogue waves |32l I33j , stock 
market crashes [3l] , optical rogue waves |35l [36] , super¬ 
fluid helium Bose-Einstein condensates [35], atmo¬ 
spheric physics [35] and opposing current flow [JD]. The 
rogue waves in plasmas may be generated due to wave- 
wave interactions [41j . and their evolution can also be 
described by the NTS equation for modulated wave pack¬ 
ets. The rogue waves were first observed in a ocean [33] . 
The space-time localized solution (rational solution) of 
the NTS equation was determined and later recognized 
as Peregrine solitons (a prototype of rogue waves). The 
latter were also experimentally observed in different non¬ 
linear dispersive media (e.g., optical fiber, plasma etc.). 
The possible mechanism to explain these waves include 
effects such as nonlinear focusing via MI in one space di¬ 
mension. For some theoretical description of rogue waves 
readers are referred to Ref 1431 and some experimental ob¬ 
servations to RefsiMi mmeiiis 

We, however, mention that the physical mechanism 
behind the formation of rogue waves is still unclear. 


however, observations indicate that they have unusually 
steep, solitary or tightly grouped profiles, which appear 
like “walls of water” mi- A large number of theoreti¬ 
cal and experimental investigations have been reported 
for the study of rogue waves in various plasma systems 
|45l IhOHHO] . Abdelsalam [58] studied the nonlinear prop¬ 
agation of solitary and freak waves in a superthermal 
plasma with positive and negative ions, ion beam and sta¬ 
tionary dust particles. More recently, Shalini and Saini 
[61] have considered the propagation of DIA rogue waves 
in superthermal plasmas and investigated the combined 
effects of superthermality and dust concentration on the 
characteristics of DIA rogue waves. 

In this work, we consider the amplitude modula¬ 
tion of lAW packets in an unmagnetized plasma whose 
constituents are inertial cold positive ions and two- 
temperature (cool and hot) electrons obeying Tsallis’ 
nonextensive q distribution. The nonlinear evolution of 
these waves are governed by a NTS equation which is used 
to study the MI of a Stokes’ wave train to a small longi¬ 
tudinal perturbation. It is shown that the wave packets 
with long-wavelength modes {k —>■ 0) are always stable, 
wheres modes with fc > I gives rise the MI of lAWs. The 
manuscript is organized as follows: In Sec. [^ we present 
the fluid model for lAWs. Using the multiple scale per¬ 
turbation technique, the expressions for different order 
harmonic modes along with the NTS equation are pre¬ 
sented in Sec. IIIIl Section [|V] describes the criteria for 
MI and Sec. [V] presents the evolution of bright and dark 
envelope solitons as well as rogons. Finally, Sec. [V^ is 
devoted to summarize and conclude the results. 


II. FLUID MODEL 


We consider the nonlinear propagation of lAWs in 
a collisionless unmagnetized plasma consisting of q- 
nonextensive distributed cool and hot electrons (to be de¬ 
noted, respectively, by the subscripts c and h) and singly 
charged inertial positive ions (with charge e and mass 
rrii). At equilibrium, the overall charge neutrality condi¬ 
tion reads nuo/niQ = 1 — /, where / = nca/riiQ is the ratio 
between the cool electron and the ion number densities. 
The dynamics of ion-acoustic perturbations is governed 
by the following set of normalized fluid equations 


duj 

~di 




( 2 ) 


d'^cj) 

dx"^ 


fric + (1 - f)nh - rii, 


(3) 


where, rij is the number density of j-species particles 
(j = c, h, i, respectively, stand for cool electrons, hot elec¬ 
trons and positive ions) normalized by their equilibrium 
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value Ujo and Ui is the ion-fluid velocity normalized by 
the ion-acoustic speed Cg = ^/ksTc/nii with ks denot¬ 
ing the Boltzmann constant. Also, (p is the electrostatic 
wave potential normalized by ksTc/e. Furthermore, the 
space and the time variables are normalized by the De- 
bye length Au = (fesTc/dTruioe^) and the ion plasma 

period Wp/ = (dTrriioe^/TOi) respectively. In Eq. ([^, 
the number densities ric and for cool and hot electrons 
(j = c, h) are given by the following g-nonextensive dis¬ 
tributions [21] 





- 1) 
ksTj 


{l/2) + l/{qj-l) 


(4) 


where Qc [qh) is a nonextensive parameter for cool (hot) 
electrons, measuring deviation from the Maxwellian dis¬ 
tribution with Qj —>■ 1. 

In the propagation of small but finite amplitude lAWs, 
one must have lep/ksTd 1. Also, as mentioned in the 
introduction that a q distribution with g < 1 is strongly 
suggested for real plasma systems or superthermal plas¬ 
mas. So, we take qc^h < 1- Thus, Eq. can be binomi- 
ally expanded (retaining terms up to to obtain from 
Eq. @ the following: 

d'^d) 

^ = l + (ci-l-di)(^-|-(c2-l-(i2)<('^ + (c3-|-d3)0^ —rii, (5) 


by taking A = A(°) + -h €^A(^^ + ■■■ = -h 

e"A("), where e <C 1 is a small positive parameter 
measuring the weakness of the wave amplitude. The wave 
amplitude is thus allowed to depend on the stretched 
(slow) coordinates of space and time as A„ = e^x 
and Tn = respectively, where n = 1, 2, 3,... (viz. 

Xi = ex, X 2 = e^x, and so forth; same for time), dis¬ 
tinguished from the (fast) carrier variables x (= Xq) and 
t (= To). All the perturbed states depend on the fast 
scales via the phase 9i = kx — cot only (w, k are the 
wave frequency and wave number respectively), whereas 
the slow scales only enter the Ith harmonic amplitude 
a[”\ Thus, A" = 'ZZ -00 in which 

the reality condition A^j^ = A|"^* is met by all the state 
variables. 

In what follows, we substitute the above expansion into 
Eqs. Q, ([^ and ([^ and equate coefflcients of different 
powers of e. Thus, equating the coefficients of e for n = 1, 
1 = 1, one obtains 

— + iku^i^ = 0, (7) 

— + ik(j)\^'^ = 0, (8) 


- (/c^-I-Cl-I-= 0. (9) 


where the coefficients Cj and dj (j = 1,2,3) are all posi¬ 
tive and given by 

ci = ^/(gc + l), C2 = i/(gc + 1)(3 - gc), 

C 3 = ^/(gc + l)(3-gc)(5-3gc), di = ^{1-f)l3{qh + l), 
d2 = ^(1 -/)/3^(g/i + 1)(3-g/i), (6) 


The solutions for the first harmonics read 

= {k'^ + Cl + ( 10 ) 

= ^(fc^ + Cl + di)(j)[^^ = (11) 

Thus, eliminating and we obtain the following 
linear dispersion relation for lAWs 

uj = k/\/ k'^ + Cl + di- (12) 


d3 = ^(1 - f)/3^{qh + 1)(3 - gh)(5 - 3qh)- 

where /3 (= Tc/Th) is the ratio of cool to hot electron 
temperatures. 


III. DERIVATION OF NLS EQUATION: 

PERTUBATIVE APPROACH 

We employ a multiple-scale perturbation technique 
[5^ IS5] to derive the evolution equation of a slowly vary¬ 
ing weakly nonlinear wave amplitude of lAW packets in a 
g nonextensive plasma. Here, we consider A as the state 
vector {A}(= m, Ui, p), describing the system’s state at 
a given position x and time t. Small deviations will be 
considered from the equilibrium state A^°) = (1,0,0)"^ 


This clearly shows the dependency of the wave frequency 
of lAWs on the nonextensive parameters gc and qt as 
well as the density and temperature ratios / and f3. It 
is seen that the wave frequency increases slowly with k, 
however, it decreases with increasing values of qj, f and 
p. We note that since ci and di are positive the linear 
stability is ensured in Eq. (12). 

For n = 2 and I = 1, we obtain a compatibility condi¬ 
tion in the form 


dTi ® dXi 


= 0 , 


(13) 


where Vg = doj/dkis the group velocity of the lAW pack¬ 
ets given by 


(A:2 


Cl 


Cl 


dl)3/2' 


(14) 
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Clearly, the group velocity of the wave becomes constant 
in the long-wavelength limit fc —>■ 0. This implies that 
in a frame moving with the group velocity Vg, the time 
derivatives of all physical quantities should vanish and 
thus one can observe a slow variation of the wave ampli¬ 


tude in the moving frame of reference. From Eq. (14) we 
also find that the group velocity decreases not only with 
the wave number k of the carrier wave, but also with in¬ 
creasing values of the plasma parameters qj {j = c,h), f 
and /3. 

Next, the expressions for the wave amplitudes corre¬ 
sponding to n = 2,1 = 1 are obtained as 


( 1 ) 


„(2) _ -2ik^ 


(15) 


M 20 ) ‘^{C 2 + d 2 )Vg - {k'^ + 3{ci + di)) 

' “ l-(Ci+di)Eg 2 ■ ^ ^ 

Finally, for n = 3, ^ = 1, we obtain equations for third 
order, first harmonic modes in which the coefficients of 
and /dXi vanish by the dispersion relation and 

the group velocity expression. In the reduced equation, 
we substitute the expressions for second order, zeroth 
harmonic modes. Thus, we obtain an equation in the 
form 


( 1 ) 


dTo 


( 1 ) 




dXo 


+ (23) 


( 2 ) . 


(16) 


Similarly, for n = 2 and 1 = 2, the evolution equations 
provide the amplitudes of the second order, second har¬ 
monic modes which are proportional to ■ The 

expressions for them are obtained as: 

(2) _ ^(22) / „,(2) _ ^(22) / ,(1)^ 2 


) , vX> = c: 


)■ 


where the coefficient of dispersion P and the nonlinearity 
Q are given by 


P= 

“ 2 5/c2 


3 fc(ci -I- di) 

2 {k? -I- Cl -I- di)®/2 ' 


Q=fc(fc2+ci+di)-3/2 


(c2+d2) (cr^+Cr^)+^(c3+d3) 


(22)) _ 


( 


k (+ 




(24) 


(25) 


4^) = 


where the coefficients are 


(d-r. 


cf = (C2 + d2) + («" + Cl + dl)C';^""^ 


( 22 ) 


The nature of these coefficients will be discussed in detail 

(17) 

in the next section. 

The total (first order) solution obtained for the electric 
potential is 

(j) ~ = (j){x, -h (j)* {x, , (26) 


^(22) _ w 
- fc L 


- (A:^ + Cl + di)2 


and 


^(22) _ (c2+d2) , (fc^ -I-Cl -I- di)^ 

03 — - 


3fc2 


2fc2 


(18) 


Now, for the zeroth harmonic modes we consider the ex¬ 
pressions for n = 2, Z = 0 and n = 3, I = 0. Thus, we 
obtain 






where the coefficients are 

= [(ci + <^i)C'3 + 2(c2 -I- ^ 2 )]; 


(19) 


( 20 ) 


^( 20 ) ^ ^ ^ ^^)2 ^ VgCi^°\ ( 21 ) 


where the variation of 4>{x, t) is assumed to be slower than 
that of exp [i{kx — ujt)] (e = 1 was formally set here, with 
the understanding that \(j)\ = |<(>4^| 1, i.e., remains 

small). The compatibility condition is now written in 
the form of a NTS equation for small and slowly varying 
wave amplitude (j){x, t) as 


d(j) 

Ih 




d(j) 

dx 


+ i "0 + QI'^P<(> = o. 


(27) 


Performing a Galilean transformation, namely ^ = x — 
Vgt, T = t, Eq. (27) can be written in its standard form 
as 




(28) 


where dispersive and nonlinear coefficients P and Q re¬ 
spectively are given by Eqs. (24) and (25). Since P and 


Q are dependent on various physical parameters viz /, 
/3, Qc and the carrier wave number k, so it is very im¬ 
portant to understand the nature of solitary structures 
from the solutions of Eq.(28) under the influence of such 
parameters. 
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IV. MODULATIONAL INSTABILITY 


We consider the a mplit ude modulation of a plane 
wave solution of Eq. (28) of the form (j) = 
where flo = is the nonlinear frequency shift 

with 00 denoting the potential of the wave pump. We 
then modulate the wave amplitude with a plane wave 
perturbation with frequency and wave number K as 
0 = (00 + 6016 *^^“*^"^ + c.c.) One thus obtains 

the following dispersion relation for the modulated lAW 
packets. 




(29) 


where Kc = -\/2|(3/P||0o| is the critical value of the wave 
number of modulation. Thus, the MI sets in for K < Kc 
and the wave will be modulated for PQ > 0, or more pre¬ 
cisely, for Q < 0 since P is always negative. In this case, 
the perturbations grow exponentially during the propa¬ 
gation of lAW packets. On the other hand, for K > Kc, 
i.e., for PQ < 0, or Q > 0, the lAWs are said to be sta¬ 
ble under the modulation. The growth rate of MI can be 
obtained from Eq. (29) as 


r = \p\K^ 


Kl 

K^ 


- I 


(30) 


Clearly, the maximum value Emax of E is achieved at 
K = Kcl\f2 and is given by E^ax = |Q||0oP- 

We note that the coefficients P, Q and hence the crit¬ 
ical value Kc as well as the growth rate E all depend not 
only on the plasma parameters qj, f and (3, but also on 
the carrier wave number k. In the range of values of k, 
for which the MI is excluded, i.e., PQ < 0 or Q > 0, 
the perturbations may develop into dark solitons which 
represent an electric potential envelope hole with finite 
or vanishing potential at the origin. Such solitons occur 
in plasmas under the effects of self-modulation. On the 
other hand, for PQ > 0, i.e., for Q < 0, the excitation 
of bright envelope solitons may be possible. Although, 
envelope solitons appear to be analytical coincidence via 
this model, and the phenomena should be clearly distin¬ 
guished [linear amplitude stability analysis in one case, 
nonlinear theory of partial differential equations (PDE) 
in the other], it has been postulated in another context 
[M] that MI is the first evolutionary stage towards the 
formation of envelope solitons. The latter in the form of 
bright type are widely observed in abundance in space 
plasmas [^I^ISS] . 

Inspecting on the expressions for P and Q it is found 
that P is always negative as ci, di > 0, however, Q can 
be positive or negative depending on the range of values 
of k and the parameters qj, f and 13. Now, for fc ^ I, 
P can be written as P ~ —Pok, i.e., P tends to vanish 
as fc —>■ 0. This is expected for an ion-acoustic mode 
where to ^ k and (f'uj/dk'^ = 0. However, Q increases as 
Q ~ qo/k in the limit k ^ 0. These are in qualitative 
agreement with earlier theoretical investigations [IHl ISS] ■ 


Also, in the limit fc —>■ 0, the real quantities po and qo 
are, in fact, both positive, given by 

" 2(ci+di)3/2’ 

1 / 1 \ 3/2 

^ 12 \Cl -I- di / [2 (c2 + ^ 2 ) - 3(ci-I-di)^] . (32) 

Thus, we find that the lAW packets with long-wavelength 
carrier wave modes are always stable under the modu¬ 
lation. In the opposite limit, i.e., for fc ^ 1, we will 
see from the numerical investigation that the waves are 
mostly modulationally unstable. 

Since P is always negative, the sign of PQ, for which 
the stability and instability of the lAW envelopes can be 
determined, solely depends on the sign of Q. So, we nu¬ 
merically investigate the nonlinear coefficient Q for dif¬ 
ferent values of the parameters qj, f and (3. Typical 
behaviors of Q with respect to k are shown in Eig. 
While the upper panel shows the properties of Q for dif¬ 
ferent values of (3 and /, the lower one shows the same 
with the variations of the nonextensive parameters qc and 
qh- It is seen that whatever be the values of the plasma 
parameters, Q is always positive in the limit fc —>■ 0 (i.e., 
at long scale perturbations of the carrier modes) as men¬ 
tioned above, whereas it becomes negative for /c > 1 as 
well as for some intermediate values of A: in 0 < A: < I, 
i.e. at small and medium scale perturbations. It is also 
evident that Q becomes maximum (hence the maximum 
growth rate of instability) in the limit A; —>■ 0 as men¬ 
tioned before except for the case where qc, qn increase 
from some positive values to the value I (see the dot¬ 
ted and dash-dotted lines in the lower panel) in which 
Q can be maximum at an intermediate value of k in the 
interval 0 < A: < I. Typical behavior of Q (in the case 
oi qc 7 ^ qh < 1) shows that initially it drops down from 
its maximum value in the domain 0 < A: ^ 0.1 and then 
starts increasing in a small interval of A; ^ 0.1, and again 
it decreases faster the larger is the wave number k of car¬ 
rier modes. However, as the values of the nonextensive 
parameters qj increase and approach a value 1 (i.e., one 
approaches the Maxwellian distribution of electrons) Q 
initially increases in a small interval 0 < k < 0.5, and 
then it decreases faster with increasing values of k. Erom 
Fig. [TJ we also observe that as the parameters /3 and / 
get enhanced, the values of Q decrease, and the ranges 
of values of k for which Q < 0 increase (See the upper 
panel). The same are true for increasing values of qc (See 
the lower panel), but decreasing values of qh [except in 
the range where A: —>■ 0 in which the value of Q increases 
with decreasing values of qh, however, the range of k for 
which Q < 0 remains unaltered (See the dotted and dash- 
dotted lines the lower panel)]. Thus, from Fig. [^we con¬ 
clude that the lAW packets with short-wavelength carrier 
modes {k > 1) are always unstable with Q < 0, P < 0 
under the amplitude modulation with plane wave pertur¬ 
bation, it can also be unstable at the medium scale with 
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k lying in 0.5 fc < 1 or more generally in 0 < fc < 1. 
In these cases the modulational instability leads to the 
formation of bright envelope solitons or highly energetic 
rogue waves (rogons). On the other hand, the lAW pack¬ 
ets can be stable (Q > 0, P < 0) for long-wavelength 
carrier modes with A: << 1 leading to the formation of 
dark envelope solitons. 




FIG. 1. Plot of the nonlinear coefficient Q of the NLS eqnation 
(281 versus the carrier wave number k is shown for different 
values of the parameters /3, /, qc and qn as in the upper and 
lower panels. Since P is always negative, the stable {PQ < 0) 
and unstable {PQ > 0) regions are corresponding to Q > 0 
and Q < 0 respectively. 


temperature of hot electrons than the cold electrons, the 
lower concentration of cold electrons than ions and more 
(less) superthermal cold (hot) electrons. 




L 


FIG. 2. The growth rate of instability F is shown against 
the wave number of modulation K with the same parameters 
P, f and qj as in Fig. The other parameter values are 
k = 0.7 and cpo = 0.5. It is seen that the lower (higher) values 
of qc, f and /3 {qn) are to suppress F with cutoffs at lower 
(higher) wave numbers of modulation. 


V. ENVELOPE SOLITONS AND ROGONS 


In what follows, we numerically investigate the growth 
rate of instability T for different values of the plasma pa¬ 
rameters Qj, f and /3 as in Fig. The behaviors of F 
are shown in Fig. The critical values Kc of the wave 
number of modulation K below which the MI sets in can 
easily be calculated. It is found that as the parameters 
/3, / and qc increase, the values of Kc increase, however, 
the same decreases with increasing values of Ph. This is 
expected as Kc ^ QjP. While the upper panel shows the 
variation with /? and /, the lower panel exhibits growth 
rates for qc and qh- We find that a small increase of the 
values of /3, / (upper panel) and qc (lower panel), leads to 
a significantly higher growth rate of instability with cut¬ 
offs at higher values of K. However, as qh increases, the 
growth rate decreases (but not significantly lower) with 
cutoff at lower K. Thus, we find that the growth rate 
of instability can be suppressed in plasmas with higher 


In the previous section we have mentioned that when 
PQ > 0, the MI of lAW packets can give rise to the for¬ 
mation of bright envelope solitons due to energy local¬ 
ization. In this case, an exact analytic (soliton) solution 
of the NLS equation (28) can be obtained by considering 


(j) = exp(i0), where ^ and 9 are real functions [see for 
details, e.g.. Refs. l63l l6^ l68l as 


4’ = 4bo sech^ 


^ 2P 


C/e- 


P^ — Ut 
Wh 


CA 
- WW 1 T 


(33) 


This solution, in fact, represents a localized pulse trav¬ 
eling at a speed U and oscillating at a frequency Hg at 
rest. The width Wh and the c onstant am plitude 4bo of 
the pulse are related to Wh = ^J2P/Q4bo■ It follows that 
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for a constant amplitude, as the parameters j3, f and 
Qc increase, the width of the envelope soliton decreases, 
however, the same increases with increasing values of Ph- 
On the other hand, for PQ < 0, the lAW packets are 
modulationally stable which may propagate in the form 
of dark-envelope solitons characterized by a depression 
of the wave potential around ^ = 0. Typical form of this 
solution of Eq. (28) is given by 


V' = fpdo tanh^ 


0 = — 
2P 


^ — Vt 
Wd 


1/2 


Vi - 2PQ^d0 I T 


(34) 


This represents a localized region of hole (void) travel¬ 
ing at a speed V. The pulse width Wd depends on the 
constant amplitude V'do as Wd = -\/2|P/Q|7V^- Clearly, 
for a constant amplitude ipdo the width of the dark en¬ 
velope solitons have the similar properties as the bright 
solitons with the variations of the plasma parameters. 
Typical profiles of the bright (at k = 0.7, 1.5) and dark 
(at k = 0.1) envelope solitons are shown in Fig. [^for the 
same parameter values as for the dotted line in the lower 
panel of Fig. It is seen that while the bright envelope 
solitons exist at small and medium scales of perturba¬ 
tions, the dark envelope solitons can exist only at the 
large scale k <^1. 

In the ranges of values of k for which the lAW packet 
becomes unstable {PQ > 0), Eq. (28) can admit highly 
energetic rogue wave solutions. These rogue waves or 
rogons, in which a significant amount of energy is con¬ 
centrated in a relatively small area in space and time, 
are generated due to the MI of the coherent lAW packets 
in the limit of infinite wave modulation period. In par¬ 
ticular, they significantly amplify the carrier wave am¬ 
plitudes, and hence increase the nonlinearity during the 
evolution of the wave packets. Now, Eq. (28) can be 
rewritten as 


■ d(l) I d^(j) 2: n 

'df + -2W ’ 


(35) 


where f = Qt and i = ^jQl2Pi. The rogue wave solu¬ 
tion of Eq. (35) that is located on a non-zero background 


and localized in both space and time can be obtained as 

iHi 


</'n(l,r) 


(- 1 )” 


Grid,?") +tgn(|,T) 
Dn{i,f) 


exp(if), 


(36) 

where G„, Pin and Dn (^ 0) are some polynomial func¬ 
tions of i and f, and n = 1,2,3, •••, denotes the or¬ 
der of the solution. The first-order rogon solution is 
obtained by considering n = 1 and it corresponds to 
the Peregrine soliton [32] in which Gi = A, Hi = 8t 
and Di = 1 + 4^^ -|- 4f^. However, superposition of 
two first-order rogue waves can lead to the generation 



FIG. 3. Plots of the bright (upper panel) and dark (lower 
panel) envelope solitons at r = 0 given by Eqs. ( |33[ ) and 
(34) respectively. It is seen that while the bright envelope 
solitons exist both at the small (e.g., k = 1.5) and medium 
(e.g., k = 0.7) scales, the dark solitons exist only at the large 
scales (e.g., k — 0.1). The parameter values of qj, f and /3 
are the same as for the dotted line in the lower panel of Fig. 

The other values are taken as 'i/'fco = V’do = 0.05, U = 0.2 
and flo = 0.5. 


of another highly energetic rogue waves with higher am¬ 
plitudes. The analytic form of these waves have been 
recently obtained by Akhmediev et al. [47] using the de¬ 
formed Darboux transformation approach. This second 
order (n = 2) rogon solution has the same form as Eq. 
(36) where the polynomials G 2 , H 2 and D 2 are given by 


G 2 = - 


I" + 7 " + 


f + 5f 2 + 


(37) 


H2 = T 


3^2 - f 2 - 2 -b 72^ 


■2 , ^ 2 f _15 


(38) 


^2 =((? + /)“ + ((? - 3 /)’ 
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64 
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-b 44f ^ -b 1 


(39) 


The profiles of these rogon solutions given by Eq. (36) 
for n = 1,2 are shown in Fig. [^ We find that for /c < 1 
{k > 1), the localization occurs at large (small) i- The 














































































first-order Peregrine soliton has been recently observed 
experimentally in plasmas [45]. However, the second- 
order rogon solution, which has been observed in water 
waves |5^, is yet to observe in plasmas. The amplifica¬ 
tion factor of the amplitude of the n-th order rational 
solution [Eq. (36)] at ^ = f = 0 is, in general, of 2n + 1. 
Hence, localized lAWs that are modeled by the higher- 
order breather solutions can also cause the formation of 
super rogue waves. In order to investigate the properties 





FIG. 4. Profiles (contour plots) of the first (n = 1) and second 
(n = 2) order rogon solutions given by Eq. (361 of the NLS 
equation (351 are shown in the upper and lower panels re¬ 
spectively. The parameter values are Qc = 0.4, = 0.3, / = 

0.2, P = 0.6 and k = 1.5. 


of these rogue wave solutions more clearly for different 
values of the parameters qj, f and /3, we plot the so¬ 
lutions at r = 0. These are shown in Fig. We find 
that while the amplitudes of both the first and second or¬ 
der rogons remain the same, their widths decrease with 
increasing values of Qc (i.e., lack of superthermal cold 
electrons), / and the temperature ratio /?. However, no 
significant change of the width is found with a small in¬ 
crease of Ph for both the rogons. 


VI. SUMMARY AND CONCLUSION 

In this work, we have investigated the amplitude mod¬ 
ulation of ion-acoustic wave (lAW) packets in an un¬ 
magnetized multi-component plasma consisting of singly 
charged positive ions and two-temperature electrons fea¬ 
turing nonextensive distributions. Using the multiple- 
scale perturbation method, a nonlinear Schrbdinger 
(NLS) equation is derived which governs the evolution 
of lAW envelopes. This equation is then used to study 
the modulational instability (MI) of a Stokes’ wave train 
to a small longitudinal perturbation. Different stable and 
unstable regions of modulation in the range of the carrier 
wave number k are found out with the variations of the 
plasma parameters, namely, the nonextensive indices qc 


FIG. 5. Properties of the first (n = 1) and second (n = 2) 
order rogon solutions at r = 0 given by Eq. ( |36[ ) of the 
NLS equation (351 are shown in the upper and lower panels 
respectively for different values of the parameters as in the 
figure. For the second order rogons an increase of = 0.6 
or P = 0.4 effectively gives almost the same results (See the 
dashed and dash-dotted lines). 


and qfi, the ratio of cool electron to ion number densities 
/ as well as the cool to hot electron temperature ratio 
p. It follows that the modulated lAW packets can prop¬ 
agate in the form of bright envelope solitons or highly 
energetic rogue waves (rogons) due to MI as well as dark 
envelope solitons where the MI is excluded. The growth 
rate of instability together with the exact analytical so¬ 
lutions of the NLS equation in the form of bright and 
dark envelope solitons as well as rogons are also obtained 
and their properties are examined numerically with the 
plasma parameters. The main results can be summarized 
as follows: 

• The fluid model governing the dynamics of small- 
amplitude lAWs in multi-component electron-ion 
plasmas is valid for the plasma parameters satisfy¬ 
ing 0 < f, P, qc, qh < 1, i-C., plasmas with more su¬ 
perthermal electrons [compared to the Maxwellian 
ones {q —>■ 1)] and lower concentration of cold elec¬ 
trons than the positive ions. 

• The group velocity dispersive coefficient P of the 
NLS equation is always negative irrespective of the 
values of k and the plasma parameters. However, 
the nonlinear coefficient Q is always negative for 
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fc ^ 1 (and any values of the parameters with 0 < 
/, 13, qc, qh < 1) and can become positive in the 
limit fc —)■ 0 as well as in a sub interval of 0 < /c < 1 
depending on the plasma parameters. Thus, the 
modulated lAW packets are always stable (unsta¬ 
ble) with long-wavelength (short-wavelength) per¬ 
turbations (carrier modes) giving rise the formation 
of dark envelope solitons (bright envelope solitons 
or highly energetic rogons). The wave envelope can 
even be unstable at a medium scale of perturbation 
with 0 < A: < 1. 

• As 13 and / increase, the value of Q decreases (hence 
decreasing the maximum growth rate of instability) 
and the range of values of k for which Q < 0 (hence 
the MI) increases. The same is true for increasing 
values of qc, but decreasing values of qh except for 
some region with fc —)■ 0. 

• As the parameters j3, f and qc increase, the values 
of the critical wave number Kc below which the 
MI sets in increase, however, the same decreases 
with increasing values of j3h- It turns out that the 
growth rate of instability can be suppressed in su- 
perthermal plasmas with higher temperature of hot 
electrons than the cool electrons, the lower concen¬ 
tration of cold electrons than ions and more (less) 
superthermal cool (hot) electrons. 

• In the range of values of k where PQ > 0, bright 
envelope solitons are shown to exist both at the 


small (fc > 1) and medium (k ^ I) scale of pertur¬ 
bations. However, dark solitons exist only at large 
scales (fc <C 1). 

• The MI of lAW envelopes can also lead to the gener¬ 
ation of rogue waves in which a significant amount 
of wave energy is concentrated in relatively a small 
area in space and time. Typical profiles of the first 
and second-order rogons are presented with respect 
to finite space and time. It is found that at an equi¬ 
librium state, while the amplitudes of both the first 
and second order rogons remain the same, their 
widths decrease with increasing values of qc (i.e., 
lack of superthermal cold electrons), / and the tem¬ 
perature ratio j3. However, no significant change of 
the width is found with a small increase of Ph for 
both the rogons. 

To conclude, the findings of the present investigation 
should be relevant both in space and astrophysical 
plasma environments where two-temperature electrons 
are the main constituents together with ions. We also 
propose that an experiment should be done to observe 
the role of higher order effects for the formation of highly 
energetic ion-acoustic rogue waves (super rogue waves). 
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